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We consider a brane scenario with a massive scalar field in the five-dimensional bulk. We study 
the scalar states that are localized on the brane, which is assumed to be de Sitter. These localized 
scalar modes are massive in general, their effective four-dimensional mass depending on the mass of 
■ the five-dimensional scalar field, on the Hubble parameter in the brane and on the coupling between 
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the brane tension and the bulk scalar field. We then introduce a purely four-dimensional approach 
' based on an effective potential for the projection of the scalar field in the brane, and discuss its 

' regime of validity. Finally, we explore the quasi-localized scalar states, which have a non-zero width 

that quantifies their probability of tunneling from the brane into the bulk. 

D 

" I. INTRODUCTION 

A lot of attention has been recently devoted to the braneworld scenario, which provides an alternative to the 
standard Kaluza-Klein compactification. In fact, this idea was initiated a long time ago by the suggestion Q that 
J> . particles could be localized on a defect embedded in a higher-dimensional space, the simplest example being the case 
0^ of a domain wall, or three-brane, in a five-dimensional bulk spacetime. It was then shown, ignoring gravity, that one 
, could localize scalars and fermions on the domain wall. 

' Recent progress was achieved on the gravitational aspect of the problem, by showing that it was possible to localize 
, massless gravitons P|, and thus to recover standard four-dimensional gravity, on a brane with tension embedded 
in an AdS bulk, with the appropriate (negative) cosmological constant A = — 6/€^, so that the brane is effectively 
\ Minkowski. In addition to the massless gravitons, one finds a continuous spectrum of massive gravitons, starting from 
• zero mass but with a weak coupling to the brane so that their contribution becomes important only on scales of the 
■ order of £ or below. 

JL The analysis of the graviton modes has also been extended to the case of a brane whose effective geometry is de 

Sitter [^Q, which is the case when the brane tension more than compensates the cosmological constant. The mode 
(~| , analysis on this background yields a massless graviton, analogous to the four-dimensional graviton, plus a continuum 
' of massive modes, which starts at {3/2)H. There is thus a gap between the zero mode and the massive modes. A 
. !^ ' recent work 5] has generalized this treatment to the case of conformally flat cosmological geometries and shown that 
the minimum gap is ^Zj^H. 

M Not surprisingly, massless bulk scalar fields have properties very similar to the gravitons. In the Randall-Sundrum 

■ - - I (RS) background, it is easy to show that, as for gravitons, there exist a localized zero mode and a continuum of 
arbitrary light states. When one allows for a small but non-vanishing mass, it turns out that the massless bound 
state disappears and turns into a quasi-localized state with finite mass and finite width, the latter quantifying the 
probability for the state to decay into the bulk . 

The purpose of the present work is to study how these properties of massive bulk scalar fields are modified when 
the brane geometry is de Sitter instead of Minkowski. Massive scalar fields are of special interest for a bulk inflaton 
model of braneworld inflation [3i H 1 where inflation in the brane is driven by a bulk scalar field, in contrast with a 
brane inflaton model, where inflation is driven by a four-dimensional field confined to the brane |lC| . 

In the present work, we study in detail the bound states of a massive bulk scalar field, localized on a de Sitter brane, 
taking also into account the possibility of a coupling between the scalar field and the tension of the brane. We show 
that, depending on the range of values for the scalar field mass M and the coupling, one can find bound states, which 
can be massless or massive. 

We then consider and extend another approach, which consists in determining an effective four-dimensional potential 
for the projection of the bulk scalar field in the brane, taking into account the five dimensional effects, via the effective 
four-dimensional Einstein equations. We are thus able to deduce an effective four-dimensional potential that combines 
the five-dimensional potential and the coupling to the brane. We then compare this approach with the mode approach 
and find surprisingly good agreement. Namely, given the five-dimensional mass, both the critical coupling allowing 
for a zero mode and that at which the bound-state ceases to exist agree well with those obtained from the mode 
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analysis. 

Finally, we also explore some aspects of the quasi-bound states (or quasi-normal modes) in the range of parameters 
where the bound states no longer exist. 

The plan of the paper is the following. In the next section, we present the model and give the equation of motion 
for the bulk scalar field. In section 3, we discuss the existence of bound states. Section 4 is devoted to an alternative 
approach based on an effective four-dimensional potential. Finally, we explore, analytically and numerically, the 
quasi-localized scalar modes, which correspond to modes escaping from the brane into the bulk. 



II. DESCRIPTION OF THE MODEL 



We consider a five-dimensional bulk scalar field (f) and a brane with a tension a coupled to this bulk scalar field. 
The corresponding action is given by 



S = J d'x^/^) (^^{R - 2A5) - ^(V0)2 - y(0)^ - 1 d^x^^)<j{,^\ 
where we have introduced a five-dimensional (negative) cosmological constant 



(2.1) 



A5 = -4. (2.2) 



In this work, we will concentrate, for simplicity, on the case of a quadratic bulk potential. 



and of a quadratic coupling, i.e. of the form 



1/(0) = iAf202^ (2.3) 



^(0) = '^o + ^0', (2.4) 



where a is a dimensionless parameter characterizing the 'strength' of the coupling to the brane. 

We will take as our background configuration, the solution of the above setup, with the usual cosmological sym- 
metries (homogeneity and isotropy along the three ordinary spatial dimensions), when the scalar field vanishes every- 
where. In this case, the bulk is effectively empty of matter and, owing to the cosmological symmetries we impose, its 
geometry corresponds to AdSs. The metric can be written in the form 

ds^ = dr^ + smh^ {r/t) [-dr^ + e^^dx^gJ , (2.5) 

where we have chosen, for simplicity, a flat slicing for the three-dimensional surfaces. 

It is possible to insert a brane in such a geometry. Assuming as usual that the bulk is Z2 symmetric, i.e. mirror 
symmetric, about the brane, the junction conditions for the metric at the brane location yield the familiar brane 
Friedmann equation 

H'-^+i^^'< (2.6) 



6 V 6 , 

In the present case, the Hubble parameter is constant, so that the brane geometry is de Sitter. The Hubble parameter 
is related to the brane position ro, which is fixed in the above coordinate system (|2.5(l . according to the expression 

= \ . (2.7) 

It is then convenient to rescale the time coordinate so that it corresponds to the cosmic time in the brane. The metric 
now reads 

ds^ = dr^ + H^e^ sinh^ {r/e)[-dt'^ + H'^e^^^d^f^j]. (2.8) 

It is also convenient to introduce, in addition to the AdS Icngthscale £, another lengthscale, defined from the brane 
tension ctq. 



4 = ( y^o ) . (2.9) 
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In the RS case, £ = £q, whereas, for a de Sitter brane, £o < £. The brane Hubble parameter is thus given by 



1 1 



Introducing the variable 

z — cosh(r/^), 

the Hubble parameter given in (|2.7I) can be reexpressed as 

1 



Hizo) 



£^z 



and the expression of £q in terms of zq is therefore 
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(2.10) 



(2.11) 



(2.12) 



(2.13) 



Let us now turn to the scalar field, which will be considered as a test field on the background configuration previously 
defined. In other words, we will ignore its backreaction on the geometry. In the bulk, the scalar field must satisfy the 
Klein-Gordon equation, which is given, for the metric (|2.8(l by 



1 



sinh^(r/£) 



1 



sinh''(r/^) 



dr (sinh'*(r/£)9^0) + M^,/, = . 



(2.14) 



In addition to this bulk equation, the scalar field must also satisfy the appropriate boundary condition due to the 
presence of the brane. It can be shown that, for Z2 symmetry, this boundary condition is given by 



1 da 

'2d6 



(2.15) 



It is clear from the two above equations (|2.14() and (|2.15() that we chose to work with quadratic expressions both for 
the bulk potential and the brane coupling in order to get linear equations for the scalar field, which greatly simplifies 
the analysis. 



III. EXISTENCE OF BOUND STATES 



The purpose of this section is to solve explicitly the system of equations H2.14l[TT^ . Since the bulk equation is 
separable, it is natural to look for solutions of the form 



(j> = u{r)ip{t), 



(3.1) 



which leads to a separation of the bulk differential equation (|2.14|l into a radial equation and a time-dependent 
equation. 



1 



sinh''(r/^) 



dr {smh'^{r/£)dr) - + 



A2 



^2 sinh^(r/£) 



u{r) = 0, 



(P d no 



^(t) = 0, 



where A2 is a separation constant. And the boundary condition H2.15j) becomes 

—uir) — ——u(r) at r = ro . 
dr £ 



(3.2) 
(3.3) 

(3.4) 



The solution of the radial equation H3.2|l can be expressed in terms of associated Legendre functions. Using the 
variable z defined in (|2.11|) . one finds 



u^{z) 



(z2 _ 1)3/4 



(3.5) 
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with 

/i2 = J-A2, V = x/M2£2 + 4, (3,6) 

The effective four-dimensional mass of the mode, which appears clearly in the time-dependent equation is thus 
given by 



m 



,2 
(4) 



In the rest of this section, we wish to investigate whether there exist bound state solutions. Bound states modes 
are defined as solutions which are normalizable, i.e., such that 



dr siVih^ r Up^{r)u* (r) — 2 / dz\/ — 1 (z) < oo. (3.8) 



Since the behaviour of the Legendre functions when z approaches 1 is given by 

Pii{z)^,^^^^^{z-l)-^'\ z^l, (3.9) 

the integral in l|3.8|l converges if is real and negative. Therefore, a bound state exists if one can find a solution H3.5|) 
with /X < that also satisfies the boundary condition (|3.4|l , which can be reexpressed in terms of the variable z as 

da 

(z2„ 1)1/2 " = ^t---o. (3.10) 
Substituting the solution H3.5|l . this gives the condition 

pM _ 3 Z „ " pAi _ 

"^-1/2 2z2-l --1/2+ (^2 _ 1)1/2^.-1/2 



({u - 2)z + a{z^ - 1)1/2) - -f - l/2)P,^3/2] = , (3.11) 



where the second expression is obtained by using the recurrence relations satisfied by the associated Legendre functions. 

The bound states are delimited by two extreme cases: the massless mode (zero mode), characterized by the value 
/i = —3/2, and the mode at the top of the gap, characterized by the mass-squared m2^j ~ {3/2)H'^, i.e., /i = 0, above 
which bound states cannot exist. For each set of values for H (or z) and M, one can thus define two critical values 
for the coupling, the lower one, a^m, corresponding to a bound state that is also a zero mode and the upper one, abs 
giving the threshold above which bound states disappear. 

Using the boundary condition H3.11|l . one immediately finds that the zero mode coupling, corresponding to /i = 
—3/2, is given by 

p-5/2 

a,^{z, M) = ~M^e . (3.12) 

^iy-l/2 

where we have used one of the recurrence relations for the Legendre functions and the definition of v in terms of 
M . Note that, for M = 0, one gets azm — 0: one recovers the usual zero mode of a massless bulk scalar field with 
no coupling to the brane. As soon as the bulk mass M is nonzero, with M"^ > 0, the zero mode coupling becomes 
negative, which means that one can find a zero mode localized on the brane only with a negative coupling. The case 
opposite to this occurs for an unstable potential, i.e., M2 < (for M2£2 > which we assume in the present paper). 

It is interesting to consider two limits concerning the Hubble parameter of the dS brane: the limit where the 
Hubble parameter is very small, H£ <^ 1, corresponding to large values of z; and, conversely, the limit where Hi 3> 1, 
corresponding to values of z very close to 1. For the first limit, one can use the asymptotic behaviour 

P,^(z)^2-^-V2 .'^^ (3.13) 

r (1 -I- 1' - ^) 

which leads to the result 

azm =i 2 - 1/ = 2 - Vm^P + 4, He^l. (3.14) 
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For the second limit, the expression (|3.9|l shows that behaves Hke the inverse of the Hubble parameter, more 
specifically 

azm^-^, H^>1. (3.15) 

Similarly, one can define the largest coupling, abs, that allows for the existence of a bound state. Using the boundary 
condition (|3.11() with = 0, it is given by the expression 

( M\ ('^-l/2)P,_3/2-(^-2)zP,_i/2 

When the mass of the bulk scalar field vanishes, Af = 0, then v = 2 and the above expression for the critical coupling 
simplifies and reads 

«b.(^,0) = ^^-^-^p5^. (3.17) 

As before, it is interesting to consider the behaviour of the maximal coupling in the two limits z ^ oo and z — > 1. In 
the first limit, one finds the same result as for ttzm, namely 



^/M2F+4, i/^<l, (3.18) 



which means that a bound state exists only for a value of the coupling very close to 2 — z^, i.e., negative for > 
(and positive otherwise). In the second limit, one finds that abs is proportional to iJ, 

abs^i^i?^, Bl~>\. (3.19) 

In Fig. n we have plotted the critical couplings azm and abs as a function of z for three bulk masses (M^^^ = 1, 
^^2^2 _ 7Vf2^2 _ This defines, for the three cases, the regions in the coupling parameter space that allow for 
a bound state. Outside these regions, there is no bound state and the four-dimensional mass of the modes has an 
imaginary part, as will be discussed later. For large z, i.e. for small H, the bound state regions become thinner and 
thinner bands around the line a = 1 — v . 

Note that, for a positive mass Af , the critical coupling abs is not always positive, but becomes negative for z > Zc(Af). 
This means that whereas for z < Zc{M), i.e., for a Hubble parameter sufficiently high, one can always find a bound 
state with a vanishing coupling, it is no longer the case for z > Zc{M) where a bound state can be found only with a 
negative coupling. 

Finally, a more general quantitative analysis is possible in the cases where Af^^^ ^ 1 and a ^ 1. Indeed, by 
linearizing the boundary condition (|3.11|l about the zero mode solution characterized by = —3/2, M — and a — Q, 
one can find the linear deviation A/i = + 3/2 as a function of the small parameters f^Pi"^ and a. Expressing A/i in 
terms of the effective four-dimensional mass of the bound state, one gets 

r2 



m^4)^/i(i/£)^+/2(7?£)^, (3.20) 



with 



and 



= '\ = VTTlPP In ^ + ^1+/ . (3.22) 

The expression (I3.2(J|I shows that, in the linearized limit, the (small) bulk mass-squared and brane coupling contribute 
additively to the effective four-dimensional mass-squared of the bound state. The coefficients in the linear combination 
depend on z, i.e. on the value of the Hubble parameter. In the limit Hi — > 0, i.e., z — > -fcx3, one finds 



(3.23) 
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a 




FIG. 1: Range of coupling values a that allow for the existence of a bound state, as a function of z (and thus of Ht) for the 
three following cases: M'^l^ = 1 (long-dashed blue lines); M^£'^ = — 1 (short-dashed red lines); M^l^ = (continuous green 
lines). In each case the upper limit corresponds to ahs{z, M) and the lower limit to a^miz, M), both values converging towards 
a = 2 - \/4 + M^P at large z. 

whereas in the opposite Umit one gets 

ml. ~ + 3^a, Ht > 1. (3.24) 

The contribution from the coupHng is thus proportional to the coupling parameter m times the largest of the two 
mass scales H and £^^. The contribution from the bulk scalar field mass is essentially the same with a very small 
variation of the coefficient. Note also that the result l|3.15|l follows immediately from (|3.24l) with m^^^ = 0. 

In the case of H£ ^ 1, it may be worth mentioning that the coupling {a/£)(fP on the brane contributes to the 
effective mass term as it is, if we rescale (/) as — > $ = by introducing an effective four-dimensional scalar field 
<I> of correct dimensions, whereas the bulk mass-squared contributes with a factor of 1/2. 

IV. EFFECTIVE POTENTIAL APPROACH 

In this section, we will use the four-dimensional projection of the five-dimensional Einstein equations onto the brane 
in order to establish an effective potential for the four-dimensional projection of the scalar field on the brane. Such 
a procedure was successfully applied in the case of a bulk scalar field with quadratic potential but without coupling 
to the brane and it was shown that the effective potential is simply half the bulk potential when H£ <C 1. The 

purpose of this section is to generalize this result to include the coupling of the scalar field to the brane. 

A. General derivation 

We will start by recalling some results derived in jl2i | following the procedure of (13i] . It was shown there in particular 
that the effective four-dimensional Einstein equations for dilaton-vacuum configurations are given by 

('^G^. = - ^'^A5^, - E,, , (4.1) 

with 



(4.2) 
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(5) A 1 



(4.3) 



where denotes the covariant differentiation with respect to the metric on the brane. Using the four-dimensional 
Bianchi identities, the covariant differentiation of Einstein's equations (|4.1|1 impUes 



(4.4) 



Speciahzing to a FLRW (Friedmann-Lemaitre-Robertson- Walker) geometry, one finds that the Friedmann equation 
on the brane, corresponding to the component (0-0) of (|4.1|l is given by 



where E 



3^ =-^ + ^ 

-Eoo- Equation (|4.4() yields 
1 



1 



1 

16^ 



E, 



(4.5) 



E + AHE 



2i^^-ir-^(.^)' + ^(.'7 



(4.6) 



where we have rewritten the terms in the brackets involving the first and second derivatives of as two linear 
combinations, one which will be easily integrated and the other one corresponding to the familiar four-dimensional 
Klein-Gordon equation. 

We are now going to assume that the brane value of the scalar field satisfies a Klein-Gordon equation that can be 
written in the form 



• + 3H(t> - 



(4.7) 



where VcS is an effective potential which we wish to determine and where J stands for a possible energy leak out of 
the brane into the bulk. Substituting this Klein-Gordon equation in Eq. H4.6|l above, one finds 



2 

E + AHE ^'^{4> + 2H^^ 4) 



1 



V'„ - -V - ■^{(j 



12 



1 

16 



J 



This strongly suggests that the effective potential, if it makes sense, is of the form 



1, 



1 



-V+^a^ - —a 



12 



16 



(4.8) 



(4.9) 



where the choice of the constant is for convenience as will be seen very soon. Remarkably, this is the same combination 
which appears in the Friedmann equation (|4.5|) . That is, if we adopt this definition for the effective potential, and 
introduce the quantity X defined by 



4 ^ 



nlX, 



the Friedmann equation ()4.5|l and the equation (|4.8|l for E reduce to the very simple system 

1 



X + AHX = J(f), 



V,s + X 



(4.10) 



(4.11) 
(4.12) 



where one recognizes the standard four-dimensional Friedmann equation with a scalar field and some extra component 
X. The second equation is a (non)-conservation equation for the extra-component X. Note that X plays the role of 
the Weyl, or dark radiation, which was identified in the simplest model of brane cosmology llj. When the energy 
outflow is zero, i.e. J = 0, one recovers the result 



X = 



(4.13) 



A non zero J means that there is an energy outflow from the brane into the bulk, which is going to feed the Weyl 
energy density. This is similar somehow to the growth of C induced by the gravitational wave emission from brane 
cosmological perturbations, which was recently analysed in [T^ . It should also be mentioned that we have not really 
proved the form H4.9|) for the effective potential. We have simply shown that this form is consistent with the system 
of effective equations in the brane. 



8 



B. Exact solutions 



It is instructive at this stage to check the above approach for known exact solutions. One of the simplest examples 
is the case of a bulk scalar field with the exponential potential [isl IT^ IT^ ITM [T9l | 



V{^) = exp 



V3 



\hi(f) 



One can find for the bulk with vanishing cosmological constant explicit static solutions, which read 



ds^ = -h{R)dT^ 



for the metric, with 



C being an arbitrary constant, and 



h{R) 



h{R) 



1 - (AV4) 



R^ - CR 



(4.14) 



(4.15) 



(4.16) 



^0 = Aln(i?). 



for the scalar field. 

A brane with a tension 



(t(0) = (To exp 



K(j) 



will undergo a cosmological evolution governed by the generalized Friedmann equation 



36 



1 - (AV4) 



R 



CR 



Since the bulk cosmological constant is zero, the effective potential is here 



1, 



1 



12 



16 



Vn 



A2 



_^ + _ 1 _ _ a;^ 



exp ( ——^\K(f> 



VcSfi exp [---j=\k(, 



It is not difficult to check that </> satisfies the effective Klein-Gordon equation ()4.7|l with 

A2 



J=-(1--)H0. 



(4.17) 



(4.18) 



(4.19) 



(4.20) 



(4.21) 



The Friedmann equation ()4.19fl can also be written in the effective form H4.11|l with the energy density X given by 



3 ( 1 - y 1 CR-'^-^ 



1- AV4 



(4.22) 



Thus, if we apply our interpretation of J that it describes energy flow from the brane into the bulk, the above analysis 
implies that the energy is actually flowing onto the brane for A^ < 2 rather than flowing out of the brane. It is 
also worth noting that, whereas the present exemple is an instructive check, it is however not very useful in practice 
because the energy exchange between the brane and the bulk is important so that the choice between the variables 
X and E is somewhat arbitrary (in fact, the expression for E in terms of the scale factor is simpler). It is really for 
the cases where the energy leak J is small that the effective potential approach makes really sense physically. 



9 



C. Quadratic case 

In the case of the quadratic bulk potential and brane coupling, H2.3l2.4fl . the effective potential suggested by the 
above analysis reads (with Eq. H2.9|) in mind), 

At the extremum cj) — the effective mass-squared is thus given by 

Ml^ = + 2^yTTW - ^, (4.24) 

where we have replaced by its expression in terms of i and H . It is interesting to note that this potential, which is 
quartic in 0, takes a double- well form for M^f^ < 0. If the present approach is valid, this implies that we may describe 
a situation of spontaneous symmetry breaking with Vcs- 

We now wish to compare the results of the mode analysis with the effective potential approach, and in order to so, 
to compare the effective mass predicted by the two analyses. In the previous section, we have identified two critical 
cases: the case when the bound state is a zero mode, i.e., m^^^ = 0, and the case when the bound state reaches the 

top of the gap, i.e., mj^^ = (9/4)77^. In the effective potential approach, the equation 
is easily solve to yield 



2z / 4z2 

,(z, M) = ^== ± + MH^. (4.26) 



Similarly, the bound state threshold is determined by solving the equation 



2 _ 2a o? AP _ 9/4 



which gives 



In both cases, we will keep only the root with the minus sign since it matches with the small coupling limit. In the 
small Hubble parameter limit, i.e., in the large z limit, the two expressions converge towards the same value 2 — as 
found previously. 

For small bulk mass and coupling, the effective potential approach can provide a reasonable approximation, even 
in the case H£ ^ 1. Indeed, neglecting the term in (|4.24|) . one gets 

Mlg ~ ^AP + 2^ Vl + H^P, (a<M£). (4.29) 

Therefore, in the limit H£ <C 1, one finds 

Mls^^M' + 2^, He^l, (4.30) 

which is exactly the same result as in H3.23|l . This strongly indicates the validity of the effective potential approach 
for Hi <C 1, at least in the case of quadratic potential and brane-coupling. It is then tempting to conjecture that this 
approach is valid for more general cases, including the case when the backreaction of the scalar field dynamics to the 
geometry is non-negligible. 

Even in the opposite limit H£ ^ 1, one finds 

Mis - + 2^ a, HI » 1, (4.31) 

which is qualitatively similar to (|3.24|) although the numerical coefficients in the linear combination are now slightly 
different. 
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V. QUASI-NORMAL MODES 



We have so far concentrated our attention on stable, bound-state modes. It is however instructive to study as 
weh the decaying, quasi-normal modes in the case where there exists no bound-state mode, although this is a more 
complicated problem. In particular, if the effective potential we derived in the previous section is valid, the term J 
which describes the possible energy leak to the bulk may be determined by studying the decay width of quasi-normal 
modes. 

Mathematically, quasi-normal modes are defined as those that satisfy the purely outgoing-wave boundary condition 
at future Cauchy horizon of AdS5. A quasi- normal mode can be found by solving Eq. (|3.11l) for complex /i with 
Re[/i] > 0, with Im[/i^] being related to the decay width of the effective four-dimensional mass m(4). Note that, if 
/i is a solution to Eq. (|3.11|l . so is its complex conjugate /i*. For definiteness, let us reserve the terminology 'quasi- 
bound-state modes' for those modes with non-zero imaginary part, and call the modes with positive real /i the 'purely 
decaying modes'. 



A. Analytical approach 



In the limit H'^i'^ <^l, one can obtain an analytic expression for the quasi-bound-state modes with least real part, 
i.e., the effective mass with smallest decay width. In this case, it is convenient to use one of the decompositions of 
the associated Legendre function into hypergeometric functions [20l |: 



1 ly fi 1 ly II 3 1 

2 ^ 2 ~ 2'2 ^ 2 ^ 2'""^ 



ly fj, v fi 
2 ^ 2'~2 ~ 2' 



1 



1 



2' 1-z^ 



with 



/tt r {—V — ii) 



-(z2-l)-(''+l)/2, ^2 = 



/7Fr(i + - m) 



(z2_l)-/2. 



(5.1) 



(5.2) 



Substituting this expression into the junction condition H3.11|l . and using z'^ — 1 ~ {H£) ^, one finds, by expanding 
in terms of H^£^, 



{ly - 2 + a)iiy - 1) + Ji(^.-2)(^.-l) + i 



r(i.)r(i-i.-^) 



(1 - i/)(i/ + 2 - a){H^fy + O {{H^£y+^) = . 



(5.3) 



We stress that this expansion makes sense only for ^^H^P ^ 1. The dominant terms are the first two terms which 
give 



2 tt2/i2 



1 



2n2 



{u-2 + a){y-l) 
V - A + a 



(5.4) 



For this equation to be consistent with the condition fi^H^i^ ^ 1, one must assume — 2 -I- a <C 1, which means that 
one must be very close to the bound state region defined before. In this case, the solutions are 



M - ±1 



- 2 



(1 - a){iy -2 + a) 



(5.5) 



In the more particular case where both M^^^ <C 1 and a ^ 1, one finds 



19 M^P + Aa 



(5.6) 



It may be worth mentioning that this result is valid irrespective of the relative magnitudes of H£, Ml and a, as long 
as they are all small compared to unity. In particular, if the quantity inside the square root is positive, the solution 
with the minus sign (i.e., /Lt < 0) is just a bound-state solution in the limit HI <C 1, given by Eq. H3.23|l . On the other 
hand, if the quantity inside the square root is negative, the dominant term of will become purely imaginary and a 
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small real part will develop in /x, as a result of the appearance of an imaginary part in fp. The imaginary part of /i^, 
however, cannot be determined from the first line of 1)5. 3|l which contains only polynomial functions of /i^. One must 
resort to the second line to compute its imaginary part, which will be given by 



r(-^)r(i + ^-^) 



(5.7) 



where one can substitute the solution (|5.5II in the right hand side. This expression agrees with U for a = 0. 
In the limit H^i^ <C M^i^ <C 1 and H'^i'^ < a < 1, one finds 



HH^lml^i^] ~ ±- (mH^ + 4a)% (5.8) 

which agrees with the Minkowski limit of in the case a = 0. 

Another case which may be studied analytically, though only qualitatively, is when a ^ 1. In this case Eq. (|3.11|) 
implies that the terms proportional to a should become OiXja). Hence, in the limit a — > cxi, /i is given by a zero of 
^1^-1/2 ('^) complex /i-plane (with Re[/x] > and Im[/i] ^ 0), which is independent of a. Thus /i converges to a 

finite value in the limit a — > c». 

Let us next discuss the purely decaying modes for which /i lies on the (positive) real axis, i.e., with no real part 
in the effective mass. They are usually subdominant in the sense that their decay widths are larger than the decay 
width of the complex quasi-bound-state modes discussed above. For Hi <C 1, the asymptotic behaviour of Legendre 
functions is given in (|3.13|) and one thus finds the zeros of the boundary condition H3.11|l when the common dominator 
r(l — [i) for all terms goes to infinity, i.e. for 

^„ = l + n, n = 0,l,2,.... (5.9) 

The other limit is ff£ ^ 1 and the Legendre functions behave according to H3.9|l . In the boundary condition H3.11|l . 
the dominant terms, when z ^ 1 are the terms proportional to P^_i/2 ^^'^ therefore the zeros of H3.11I) are given by 

Mn = ^ + V4 + M2 + n = 0,l,2... (5.10) 

where n is a positive integer, corresponding to cases where the Gamma function in the denominator of Pl^_i/2 blows 
up. These zeros thus now depend on M but not on the coupling a or on z (provided however z ^ I). 



B. Numerical approach 



In order to study the quasi-normal modes in more details, we have solved numerically the following equation for 
the complex number fj,, 

F{ti- M, a, z) = -(p + 1/2)-^^ + 2)z -f a{z^ - if' ^ , (5.11) 

for various values of the parameters A/, a and z. As mentioned earlier, if ^ is a solution, so is /i*. Hence, without loss 
of generality, we may confine our search for /i on the first quadrant (i.e., Re[/x] > and Im[/i] > 0) on the complex 
/j,-plane. 

Let us start our numerical exploration with the case of a vanishing coupling a = 0. In Fig. |21 we have plotted the 
numerical solutions for /i in the complex plane for various values of the bulk mass M . As can be seen in the figure, 
we have obtained several branches of solutions which evolve continously as Mi varies. As Ml increases, the modes 
on these branches migrate away from the real axis. The branch closest to the vertical axis is the most important 
dynamically because its imaginary part is small. 

This branch is connected to the origin of the complex plane in the limit Ml —t and the corresponding quasi- 
normal modes have been computed analytically just above for Ht <^ 1 and M'^t'^ ^1. As a check, we can compare 
the analytical results given above with our numerical solutions. As noted before, we perform our search in the first 
quadrant of the complex /x-plane, which means Im[m^^^] < 0. First, we have checked that the real part of Jti^^^ is 

indeed very close to the analytical value M^/2. For the imaginary part, we have plotted in Fig. 13— Im[m^^-|]/(M^i?^) 
as a function of Mi in two cases, z = 50 and z = 1000. In both cases, starting from high values of M£, one approaches 
the analytical value of 7r/16 as M£ decreases, as suggested by (|5.8|l . But, whereas in the case z = 1000, the dots 
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FIG. 2: Evolution in the complex /^-plane of the quasi-normal modes when the mass M varies (with a = 0, z — 10). The 
highest points for each branch correspond to A'I£ — 10. The increment between two adjacent points is A{M£) = 1. 
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-Im[m'^ 



(4) 



]/(M*£'^) as a function of M£ {z — 50 for the lower points, 



1000 for the upper points). The curves 



correspond to the analytical estimate 15.71 



continue to approach 7r/16, there is a sudden change in the evolution for z = 50. This is simply due to the fact that, 
even if Hi ^ 1, the small values of M£ become of the same order of magnitude as Hi and (|5.8|l is no longer a good 
approximation: one needs the more general expression H5.7|l . The two curves in the figure correspond to the analytical 
estimate given by (|5.7|l and one sees that the numerical values converge towards these curves for low masses. 

Our numerical treatment allows us to continue this quasi-normal branch beyond the regime of validity of the 
analytical calculations, and in particular for very high masses. Moreover, as mentioned before and as illustrated in 
Fig. 121 we have found other solutions for /i, away from the real axis, which correspond to additional branches. Each 
of this branch starts from the real axis after a critical mass threshold has been reached and then evolves away from 
the real axis as Mi increases. 

Let us now consider the situation when the coupling is allowed to vary. We have plotted in Fig. ^ the evolution of 
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FIG. 4: Evolution in the complex /i-plane of the quasi-normal modes when the coupling varies, in the clockwise direction, from 
to 7 (with Ml = 0, 2 = 10). The increment between two points is Aa = 0.2. 
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FIG. 5: Variation in the complex /i-plane of the quasi-normal modes when the coupling changes (for M£ = 5, z = 10). The 
increment between adjacent points is Aa — 1. The three big dots (in black) correspond to the modes for a — 0. The points on 
the left of each black dot (in red) correspond to negative couplings and those on the right (in blue) to positive couplings. 



the quasi-normal mode with the increasing (positive) coupUng in the case M£ — 0. What can be observed is that the 
mode, after going away from the real axis, tends to come back towards it. 

In Fig. 13 we have also plotted the same evolution, i.e. with increasing coupling, in the case M£ — 5 where several 
quasi-normal modes are already known to exist in the complex /i-planc for vanishing coupling. We have no analytic 
expression for this case, but one can observe each of the modes follows a curve analogous to the M£ = case. We 
have also added the evolution corresponding to increasing negative coupling. For the first branch, one finds that 
the mode converges towards the origin of the complex plane, which could be expected since at some critical value of 
the coupling one will enter into the localized mode region. For the second branch, one observes that the modes at 
large negative couplings converge towards the modes of the previous branch at large positive couplings. The same 
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behaviour occurs for the third branch. 



VI. CONCLUSIONS 



In the present work, we have studied the modes of a massive bulk scalar field that are localized or (quasi-localized) 
on a brane with a de Sitter geometry, characterized by a Hubble parameter H . We have allowed for a coupling 
between the brane tension and the bulk scalar field, quantified by a dimensionless parameter a. 

Although the bound state of a massless bulk scalar field without coupling to the de Sitter brane is a zero mode, 
i.e. its four-dimensional effective mass TO(4) vanishes, this is no longer the case for a non-zero bulk mass M or a non 
zero coupling. One thus finds in general a massive bound state, whose four-dimensional mass depends on M, a and 
H. We have computed explicitly this dependence in limit of a small bulk mass and small coupling. 

The mass-squared m^^^ of the bound state is always comprised between zero and {3/2)H^, value above which one 
finds a continuum of modes. We have obtained explicitly, for any value of M and H, the lower and upper critical 
values for the coupling corresponding to this range of possible zero modes. 

Another approach consists in trying to find a purely four-dimensional description of the scalar field including the 
effect of the bulk. We have generalized previous such procedures to take into account the coupling, and defined an 
effective potential for the value on the brane of the scalar field. For small M and a, the second derivative of the 
effective potential yields an effective mass-squared which is in excellent agreement with the rigorous mode approach 
for small H and in qualitative agreement for large H. We have also demonstrated that a class of known exact solutions 
with a bulk scalar field indeed conforms to the effective potential approach. 

Finally, when the three parameters M, a and H are not in the region allowing for bound states, one finds modes 
whose effective mass has an imaginary part. This means that the corresponding states cannot stay localized on the 
brane but will escape into the bulk. We have explored numerically the dependence of the real and imaginary part 
of the four-dimensional mass on the values of the bulk mass and of the coupling and we have given an analytical 
estimate of the quasi- localized modes for small values of the Hubble parameter, i.e. Hi ^ 1. 



Acknowledgments 

We would like to thank K. Koyama for useful communications, in particular, for bringing our attention to exact 
solutions that conform to the effective potential approach. This work is supported in part by Monbukagaku-sho 
Grant-in- Aid for Scientific Research (S), No. 14102004. This work was initiated when MS was visiting the gravitation 
and cosmology group (GReCO) at lAP, Paris. MS is grateful to the Univcrsite Paris 7-APC for financial support and 
to the GReCO members for warm hospitality. 



[1] 


V. 


ar] 


[2] 


L. 


[3] 


J. ( 


[4] 


D. 


[5] 


A. 


[6] 


S. ] 


[7] 


Y. 


[8] 


N. 


[9] 


Y. 


[10] 


R. 


[11] 


P. 1 


[12] 


K. 


[13] 


T. 


[14] 


D. 


[15] 


R. 


[16] 


H. 


[17] 


D. 


[18] 


D. 


[19] 


K. 


[20] 


see 



A. Rubakov and M. E. Shaposhnikov, Phys. Lett. B 125, 136 (1983); K. Akama, Lect. Notes Phys. 176, 267 (1982) 
Xiv:hep-th/0001113 . 

Randall and R. Sundrum, Phys. Rev. Lett. 83, 4690 (1999) arXiv:hep-th/9906064 . 

Garriga and M. Sasaki, Phys. Rev. D 62, 043523 (2000) arXiv:hep-th/9912118 . 

Langlois, R. Maartens and D. Wands, Phys. Lett. B 489, 259 (2000) arXiv:hep-th/ 0006007| . 
V. Frolov and L. Kofman, arXiv:hep-th/0209133 



L. Dubovsky, V. A. Rubakov and P. G. Tinyakov, Phys. Rev. D 62, 105011 (2000) arXiv :hep-th /0006046| . 
Himemoto, T. Tanaka and M. Sasaki, Phys. Rev. D 65, 104020 (2002) arXiv:gr-qc/0112d27'r 
Sago, Y. Himemoto and M. Sasaki, Phys. Rev. D 65, 024014 (2002) arXiv:gr-qc/01 04033. . 

Himemoto and M. Sasaki, Phys. Rev. D 63, 044015 (2001) arXiv:gr-qc/0010035 . ' 

Maartens, D. Wands, B. A. Bassett and I. Heard, Phys. Rev. D 62, 041301 (2000) arXiv:hep-ph/9912464l. 
Binetruy, C. Deffayet, U. Ellwanger and D. Langlois, Phys. Lett. B 477, 285 (2000) arXiv:hep-t h/9910219, . 
i. Maeda and D. Wands, Phys. Rev. D 62, 124009 (2000) arXiv:hep-th/0008188 . 

Shiromizu, K. i. Maeda and M. Sasaki, Phys. Rev. D 62, 024012 (2000) arXiv:gr-qc/9910076l^ 

Langlois, L. Sorbo and M. Rodriguez-Martinez, Phys. Rev. Lett. 89, 171301 (2002) arXiv:h ep-th /0206146| . 
G. Cai, J. Y. Ji and K. S. Soh, Phys. Rev. D 57, 6547 (1998) arXiv:gr-qc/9708063 . 

A. Chamblin and H. S. Reall, Nucl. Phys. B 562, 133 (1999) arXiv:hep-th/9903225(; 

Langlois and M. Rodriguez-Martinez, Phys. Rev. D 64, 123507 (2001) arXiv:hep-th/01 06245| . 
Langlois, to be published in Prog. Theor. Phys. Suppl. No. 148 arXiv:hep-th/0209261 . 
Koyama and K. Takahashi, arXiv:hep-th/0301165 

e.g., §3.2 (21) of H. Bateman, "Higher Transcendental Functions, Vol. I" (McGraw-Hill, New York, 1953). 



